Noise Suppression for Micromechanical Resonator via Intrinsic Dynamic Feedback 



o 
o 

(N 
> 

o 
in 



i 

c 

o 
o 



(N 
> 

O 

p 

o 
l> 
o 



X 



H. Ian0 Z. R. Gong, and C. P. Suifl 

Institute of Theoretical Physics, The Chinese Academy of Sciences, Beijing, 100080, China 

(Dated: February 1, 2008) 

We study a dynamic mechanism to passively suppress the thermal noise of a micromechanical 
resonator through an intrinsic self-feedback that is genuinely non-Markovian. We use two coupled 
resonators, one as the target resonator and the other as an ancillary resonator, to illustrate the 
mechanism and its noise reduction effect. The intrinsic feedback is realized through the dynamics of 
coupling between the two resonators: the motions of the target resonator and the ancillary resonator 
mutually influence each other in a cyclic fashion. Specifically, the states that the target resonator 
has attained earlier will affect the state it attains later due to the presence of the ancillary resonator. 
We show that the feedback mechanism will bring forth the effect of noise suppression in the spectrum 
of displacement, but not in the spectrum of momentum. 

PACS numbers: 85.85.+j, 85.25.Cp, 45.80.+r 



I. INTRODUCTION 



Recently, interests have been generated on cooling 
techniques for mechanical systems at nano- and micron- 
scales pj. Among them, the typically employed is the 
feedback cooling technique where an external feedback 
circuit is responsible for detecting the motion of the tar- 
get and feeding a counteracting force against this motion; 
through a general decrease of magnitude in the density 
noise spectrum, it was shown that the feedback can ef- 
fectively reduce the fluctuation of the target and provide 
a cooling mechanism 0, 0| • Some experiments based on 
the models that contain the feedback loops have been car- 
ried out in the past few years. A few are directed towards 
the cooling of micron- to nanometer-size mechanical res- 
onators, aiming to reach a macroscopic quantum mechan- 
ical ground state and serving as a powerful manifestation 
of quantum mechanical effects @,H, 0, Hi- Other experi- 
ments succeed in slowing down the motion of micron-size 
mirrors through the radiation pressure of an optical field 
in a Fabry-Perot (FP) cavity, aiming to reach a noise 
level and equivalently an effective temperature that are 
pertinent to the employment of high-precision detection 
of gravity waves @, [Hj, [HI • 

The forementioned implementations of feedback cool- 
ing through reduction of noise fluctuations invariably rely 
on an electrical circuitry external to the target system to 
be cooled. The controller here is usually fixed and at- 
tracts or repels the resonator through either electrostatic 
Coulomb force or Lorentzian force. If such an external 
detection-control unit could be eliminated in favor of a 
mechanism with self-detection of and self-adjustment to 
the target's thermal motion, we call the mechanism "self- 
cooling" [l3|, [l5| . Devices implementing this self-cooling 
use less components and are free from the reliance on an 



'Electronic address: hian@itp.ac.cn 
t Electronic address: 
URL: http://www.itp.ac.cn/suncp 



suncp@itp.ac.cn 



external circuitry and hence prone to less noise sources. 

In one case [13j . an augmented cavity along with an 
extra optical cavity field is established on the other side 
of the mirror, in addition to the regular FP cavity, so 
as to counteract the radiation pressure from the original 
cavity field. This extra field cushions the motion of the 
pressure mirror and plays the role of feedback. In an- 
other case fl3.[l5j. a set of Josephson junctions behaving 
as a qubit, serially connected to a mechanical resonat- 
ing beam, serves delayed supercurrent into the circuit 
according to the magnetic flux through the circuit loop. 
The magnetic flux is controlled by the vibrating motion 
of the beam, which in turn is controlled by the magnetic 
field generated by the current feed. Such a mutual de- 
pendence furnishes a self-feedback mechanism. It should 
be pointed out that both of the self-feedback setups re- 
quire delayed feedback, which assumed a priori a non- 
Markovian approximation that explicitly depends on the 
history of the target's motion. In these phenomenolog- 
ical treatments, the cooling target either couples itself 
to a static controller and makes itself prone to the noise 
stemmed from the feedback, or couples to a mechanically 
static detection construct and receives manually delayed 
feedback. 

Hereby, we present a dynamic model based on an in- 
trinsic mechanism with non-Markovian feedback, which 
is obviously free from an external feedback loop and does 
not rely on a presupposition of historical dependence. 
This mechanism is illustrated by a simple mechanical sys- 
tem in which the target is modeled by a harmonic oscilla- 
tor and attached to a dynamic controller, which is a rela- 
tively heavier resonator, through a spring. The target is 
controlled by an intrinsic feedback through the dynamics 
of coupling: earlier positions and velocities of the target 
affects the motion of the controller and this influence is 
subsequently fed back to the target. Consequently, the 
accumulation of earlier states of the target will affect the 
state of itself later. With proper parameter setup, the 
target essentially experiences a resistance and decelerates 
its motion; its displacement variance is shrunk, noise sup- 
pressed and effective temperature cooled down. The lack 
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FIG. 1: (Color online) The diagrammatic figure shows the 
arrangement of the three springs and the two masses. From 
left to right, they are: the spring of constant k, the target 
mass m, the spring of constant g, the ancillary mass M, and 
the spring of constant K. G(t) is the harmonic driving force. 

of a specific detection device for the motion of the target 
resonator and an external feedback circuit characterizes 
the intrinsic nature of the mechanical feedback. Our nu- 
merical analysis shows the existence of a noise suppres- 
sion capability of our scheme, e.g. the variance of dis- 
placement can be reduced to 0.04 x 10 _21 m 2 , and there- 
from a cooling capability under a practical setting acces- 
sible in current experiments. We note that the scheme is 
theoretically illustrative through its simple model setup 
yet widely applicable because the general oscillator sys- 
tems can be extended to quantum bosonic systems and 
other cases. In fact, a similar model and mechanism has 
been proposed to actively cool down the torsional vibra- 
tion of a nanomechanical resonator through spin-orbit 
interactions 

The model will be explained in SecHIl and its delay 
function then derived a posteriori to examine its non- 
Markovian dependence. The complete solution of the 
system dynamics is given in Sec lIII| with which we will 
derive the density noise spectrum and calculate the theo- 
retical noise suppression rate. The associated numerical 
results will be presented in Sec lIVl given various param- 
eter setups. The analysis is extended to the domain of 
momentum noise in SecfVl 



II. INTRINSIC FEEDBACK BY COUPLING 
DYNAMICS 

A. The Model 

Our model setup (see FigfT]) comprises two masses and 
three springs. The two masses are denoted by to and 
M, respectively. The mass to is the target and typically 
lighter whereas the mass M serves as an ancillary con- 
troller and is relatively heavier. The three springs are 
denoted by their Hooke's constants k, g and K, respec- 
tively. The spring of constant k attaches the lighter mass 
to to the fixed wall on the left and the spring of constant 
K attaches the heavier mass M to the fixed wall on the 
right. The spring of constant g strings the two masses to- 
gether. Such a setup, intuitively, grants the heavier mass 



M the function of a suspension system and a medium 
for the feedback. The symbol G{t) represents an exter- 
nal driving force which is necessary for the discussion of 
cooling but can be deemed zero for the present. 

The connected springs will give rise to mechanical vi- 
brations of the masses. We let Q = wjje + gj/m de- 
note the effective mechanical resonance frequency for 
the mass to, assuming the other mass M is fixed, and 
Cl = \J (K + g) JM the equivalent for the mass M, as- 
suming the mass to is fixed. Besides these mechanical 
vibrations, we assume each of the masses experience a 
frictional damping and we let 7 denote the damping co- 
efficient for the mass to and T that for mass M. 

Then according to the setup above, the coordinates 
of the two masses obey a coupled system of classical 
Langevin equations 

x + jx + cj 2 x — —Q = f (1) 

TO 

Q + TQ + Ci 2 Q- j-x = F (2) 

where x is the coordinate of the mass to and Q that of the 
mass M. f and F on the right hand side of the equations 
denote the random thermal noise generated by the mass 
to and M, respectively, due to their frictional damping. 

The frictional damping terms as dissipation and the 
random noise terms as fluctuations constitute the total 
Brownian thermal force in the classical Langevin formal- 
ism. This thermal force induces an thermal environment, 
the effect of which is divided among the two terms f(t) 
and F(t) according to the fluctuation-dissipation rela- 
tions [r| 

(/(t)/(i')> = 2k B T^S(t-t'), (3) 

TO 

(F(t)F(t')) = 2k B T^S(t-t'). (4) 

For now we assume the two masses are independently in- 
teracting with two thermal environments. That is, there 
are no correlations between the fluctuations of the two 
masses 

Actually, we can realize from Eq.(fT]) and Eq.([3]) that 
the motion of the mass to is resisted by a frictional force 
7, which in turn is transduced into thermal energy and 
heats up itself. This process, however, is mediated by 
the mass M that stands between to and the fixed wall 
through the term (g/M)Q. If there were not the mass 
M, the kinetic oscillation of to would be instantly reacted 
by the surrounding springs of Hooke's constant k and K. 
With the presence of the mass M and the extra spring 
of constant g, the oscillation of the mass to will first 
squeeze the spring of constant g, and then the squeezed 
spring will release and push the mass M to the right. It 
follows that the spring of constant K will be squeezed 
successively. The mediating mass M breaks the original 
single spring into two and permits these two springs to 
stay in different states, squeezed and released, and hence 
essentially delays the reaction of the springs. This cas- 
caded process is then reflected by the wall and executed 
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FIG. 2: Schematic illustration of the feedback mechanism 
shows the controller M receives signal x(t) and outputs signal 
Q(t) according to the motion of the target m. 



in a reversed order to the left; the delayed reaction of the 
springs of constant g and K acts back onto the oscillating 
mass to. The delayed reaction depends on the oscillating 
velocity of the target mass m and can thus be considered 
a feedback onto itself. 

The entire process can be regarded as a feedback loop 
from the view of control theory through the flow diagram 
shown in FigO The target resonator m is the system to 
be controlled and the ancillary resonator M becomes the 
controller which detects the signal x(t) as its input and 
feeds the signal Q(t) as its output. The cycle time of the 
loop corresponds to the delay of the controller-to-system 
reaction. The dependence of this delayed reaction on 
the tunable parameters, mainly those spring constants, 
allows us to control this self-feedback precisely to coun- 
teract the oscillating motion of the target mass. 

The weakened motion is converted to a reduction of ef- 
fective temperature of the target mass through an equiv- 
alent relation between the autocorrelation of the tar- 
get's displacement and temperature, derived from the 
fluctuation-dissipation relation Eq.([3]). Under normal 
circumstances, that is, when the coupling mass M were 
not present, this equivalent relation can be written as 



(x(t)x(t r )) = 2Trk Bl T£ {t - t') 



(5) 



where 7 is defined as in Eq. ([T]) and £0 is a function to be 
given explicitly. Experimentally, the statistical variance 
of the coordinate is measured and the above relation is 
used to compute an effective temperature. 

Increasing the damping coefficient 7 will certainly in- 
crease the enveloping rate of the oscillating motion of 
x(t). However, this reduction of motion cannot lead to 
an equivalent noise suppression effect upon the target 
mass to because of the constraint imposed by the friction- 
induced fluctuation relation described above. In fact, ac- 
cording to this relation, increasing damping rate leads to 
more thermal dissipation of the system. That is to say, 
simply enlarging frictional force merely results in a heat- 
ing effect upon the target system. Our aim, therefore, 
is to reduce the temperature of the target mass by re- 
ducing its coordinate variance through the feedback that 
does not simultaneously increase the damping coefficient. 
After we derive the explicit feedback response below, we 
will show the feedback is actually shrinking the variance 
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B. The Feedback 

Generally, a feedback external to an oscillating system 
has the effect of adding an extra driving force term on 
the right hand side of the equation of motion, i.e. 



y + jy + uj y = f + 



(6) 



where we have let y denote the coordinate of a general dy- 
namic system, Tro the feedback force, and 7, u>, / terms 
of similar meanings to those defined in Eq. |T]) . If the tar- 
get were to be cooled, the motion of the system should 
be slowed down. In other words, behaving as a func- 
tion of the target's velocity, the force should have an 
equivalent effect of increasing the damping coefficient, 
but do not increase the fluctuations. Besides, the feed- 
back force should be dependent on the target's position 
and the noise source. 

Summing up these requirements, we expect J-fb to be 
a function of y(t), y(t) and /. To overcome the generic 
fluctuation-dissipation relations and hence accomplish an 
efficient noise suppression, we shall use a non-Markovian 
type feedback: Fr, not only depends on the current value 
of the velocity and the position of the target, but also 
their historical values at past times. If we represent the 
historical dependence by a time-derivative &K(y, y, r) /dr 
and weigh the contribution of the histories by a delay 
function h(t — r), where t stands for the current time 
and t for the time in the past, the feedback force can be 
written as an integral with respect to r, 



dT dK(yy,r) 
dr 



(7) 



The above formula shows the mathematical character 
of a general feedback force. Inversely, any function that 
can express the same character should be considered a 
feedback force. Therefore, we can verify the dynamic 
response of the coupling mass M in our model as an ef- 
fective feedback mechanism by finding the corresponding 
specific expression for K{y, y, t) and hit—r). To do so, we 
solve Eq.([2]) by using Fourier transforms and integration 
by parts (the detailed derivation is given in Appendix [A)) 
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denotes the inhomogeneous part of the equation, i.e. the 
external driving force to M. In the solution, 



m = "i exp 
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is the delay function that we look for and 
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K(x,t) = 



V4ft 2 - r 2 



(x,t) 



(12) 



correspondingly shows the historical dependence. 

If we plug the solution into Eq.fl]), i.e. reduce the 
degree of freedom of the variable Q(t) in the equation 
by substituting with its formal solution, we arrive at an 
integro-differential equation of only a single variable x(t), 
which is comparable to the feedback-containing equation 
of motion that appears in previous literature [la fl5| . 
However, in the latter case the delay function is assume 
to be a non-Markovian approximation 
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(13) 



and K(y,y,r) is left unknown. Here, we have explicitly 
implemented a self-feedback mechanism that exerts force 
back onto the mass m after certain delay through the 
use of the coupling mass M. The sinusoidal factor in 
the delay Eq. lfTTj) illustrates the damped oscillating mo- 
tion of the mass M. The implicit time derivative of x(t) 
within the integrand implies up to an equivalent effect an 
additional friction to the motion of the mass m, which 
damps the oscillation without increasing the target's fluc- 
tuation. The non-integral term is Markovian and has the 
same order as lj 2 x. Though it does not appear in Eq.([7]), 
this term effectively reduces the oscillating frequency and 
shall not counteract the feedback effect. Therefore, both 
terms impose noise suppression effect to our target. The 
only limitation, however, is the thermal fluctuation F(t) 
from the coupling mass M itself and it will result in a 
limit for the suppression because of the competition be- 
tween this fluctuation and the effective feedback. 



III. EXACT SOLUTIONS OF THE LANGEVIN 
EQUATIONS 

In order to examine the validity of the above proposed 
mechanism and verify the efficacy of the cooling rate, we 
find the analytic solutions of the Langevin equations ((T|) 
and ([2]) to reflect the displacement of the target m as a 
response to its own thermal noise, the motion of the an- 
cillary mass M and the thermal noise of mass M. From 
the response function, we shall derive the effective damp- 
ing coefficient and vibrating frequency of the target mass 
m as well as the autocorrelation function of its coordi- 
nate. The noise spectrum and the total noise fluctuation 
are then defined upon these derived quantities. 



A. The Noise Spectrum 

The displacement spectrum can be written as the sum 
of responses of the noise terms (the derivation is given in 
Appendix [Bj) 



x(u) = Lf(u)f(u) + L f (uj)F(uj). 
Note that the two different susceptibilities 
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UJ 2 + iu)J e 



and 



L f (uj) 



g/m 



[Cl 2 — to 2 + iuT] [uj 2 — uj 2 + iujj e 



(14) 
(15) 
(16) 



reflect the system's different responses to thermal fluctu- 
ations from the mass m and the mass M. Here, we have 
defined the effective vibrating frequency of the mass m 
to be 
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(17) 



mM (n 2 - uj 2 ) 2 + uj 2 T 2 
and the effective damping coefficient of the mass m to be 



7eM = 7 + 



mM (n 2 - lo 2 ) 2 + lu 2 T 2 



(18) 



Eq. (Tl"4"|) means that the mechanical susceptibility or re- 
sponse function of the target m is adjusted because of the 
dynamic coupling of the target mass m to the ancillary 
mass M. The first term represents the direct effect of the 
thermal bath acting on the target mass m; whereas the 
second term represents the indirect effect of the thermal 
bath onto the mass m through the mediating mass M 
and the coupling between the two masses. 

The positivity of the second term in Eq. (Tr8")) has as- 
serted our expectation of increasing the damping rate 
without increasing thermal force. We shall also note 
that an additional noise source F(u>) is imposed onto 
the adjusted susceptibility Eq. (fT5)l due to the dynamic 
coupling. But seeing that it is divided by a frequency- 
squared term, we expect it to be negligible when the tar- 
get mass m is not resonating at a frequency close to that 
of the mass M. 

The density noise spectrum (DNS), which is defined 
by the equation 



S x {u>) 
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dto' (x(lu)x(u)')} , 



can thus be computed from Eq. (fl4j) 
2k B T le 



S x (uj) 



m 



w 2 7 2 ' 



(19) 



(20) 



Comparing the above equation with the case when M is 
absent, i.e. when 7 e is degenerated to 7, we observe a 
general suppression at the noise peak and spreading of 
the noise spectrum. 
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B. Noise Suppression 

The observable effect of the noise sources on the mo- 
tion of the mass m is equivalent to the variance of the 
displacement of the mass m in time domain, which is de- 
fined as the average of the entire noise spectrum, i.e. the 
integral of the DNS of x, 



duj S x (uj) 

-co 

7 rji p OO 



mir 



(uj 2 — uj 2 ) 2 + uj 2 y t 



(21) 



;dLU. (22) 



The integral is computable after we approximate the ef- 
fective damping rate and vibrating frequency by truncat- 
ing their expansions shown in Appendix [C] As a result, 
the variance is a variable of temperature and the spring 
constants 
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(23) 



The noise fluctuation depends on the three spring con- 
stants k, K and g as its parameters, and is independent 
of the masses m and M and the damping rate 7 and V of 
the resonators. We shall notice that the noise suppres- 
sion effect for the target resonator is always present for 
all values the spring constant g takes. This ideal result is 
due to the non-Markovian feedback we derived in SeclIIl 
which always increases the effective damping while re- 
taining the same fluctuations. 

Eq. ((23|) shows a complex relation between itself and 
its three parameters of spring constants. To illustrate its 
behavior, we focus on its relation with the spring constant 
g. It is a monotonic decreasing function of g and its two 
limiting values are 



<* 2 > 
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\ / I g — ^00 



k B T 
k 

k B T 
k + K' 



(24) 
(25) 



which coincides with our expectation that enlarging the 
constant g will render the feedback more effective due to 
the enlarged feedback amplitude Eq.(|TD]), 

The symmetry between the target mass and the ancil- 
lary mass in the model setup allows us to compute the 
variance of the ancillary mass following the same method- 
ology 



• + k 



K) + kK 



(26) 



The limiting value of (Q 2 ) at g — > 00 is identical to that 
of (x 2 ), which shows that under the extremal case where 
the two oscillating masses combines into one by a rigid 
body the limiting behaviors of the two bodies become the 



IV. NUMERICAL ANALYSIS OF NOISE 
SUPPRESSION 

To show completely the noise suppression behavior of 
the two oscillating masses, the numerical analysis is sep- 
arated into two parts with each part for each extremal 
end of the values of the spring constants whereas the 
other parameters are set to laboratory accessible values 
for common micromechanical resonators. The first case, 
the identical case, is where the springs attaching the two 
masses to the walls share the same Hooke's constant, i.e. 
k = K. We examine how the noise suppression behav- 
ior is affected by varying the value of the constant g of 
the middle spring. The second case, the large detuning 
case, occurs when the springs attaching the two masses 
to the walls take vastly different values of their Hooke's 
constant. Again we examine the noise suppression limit 
for different values taken for the constant g of the middle 
spring. 



A. The Identical Case 

We assume the springs with one end fixed to walls have 
k = K = lNm~ x . The lighter target mass has m = 
1 x l§~ s kg and its factional damping rate 7 = 0.1s . 
The heavier ancillary mass has M = 1 x 10~ 6 kg and 
its frictional damping rate T = 4s" 1 . When free from 
the stringing spring g, the target mass will oscillate at 
a natural frequency of 10kHz and the ancillary mass at 
1kHz. The system's initial temperature is set to room 
temperature T — 295K. 

We first look at the density noise spectrum of the ancil- 
lary resonator when the middle spring is set to have its 
constant g = O.OliVm -1 , O.lJVm -1 , lNm' 1 , lONm" 1 
and lOOA^m -1 as shown in Figf3] The noise peaks at the 
frequencies uj c = 1005Hz, 1044Hz, 1223Hz, 1376Hz and 
1404Hz. 

We notice that when tuning the spring constant g, 
not only the peaking frequency is shifted to the right, 
the peak amplitude is reduced along with the increased 
value of g. This proves a general suppression in noise 
and an equivalent cooling effect to the system. The de- 
tails of the spectrum and the spread width can be shown 
more apparently when we cluster the peaks together with 
a common frequency lu c for their corresponding peaking 
frequencies, which is shown in Fig[?J The total noise fluc- 
tuation reached after suppression can be computed from 
the area under each curve in the figure, using Eq. (|22[) . 
The attenuated noise levels corresponding to the 5 val- 
ues of the spring g are, respectively, 4.03, 3.73, 2.72, 2.12 
and 2.04 times a common factor of 10 _21 m 2 . 

The cooling effect for the target mass m is more obvi- 
ous if we examine the clustered peak plot of the target 
mass' density noise spectrum shown in FigO However, 
differing from the behavior of the ancillary mass shown in 
FigEl the target mass is resonant at two peak frequencies 
for each of the 5 values of the spring constant g. Among 
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FIG. 3: (Color online) Plot of DNS of the ancillary resonator 



M with parameters: k ■■ 
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top to bottom are plotted from g 
liVm -1 , 10iVm -1 , lOOiVrn -1 . 
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FIG. 4: (Color online) Rescaled plot of DNS of the an- 
cillary resonator M with parameters: k = K = lNm~ , 
M = 1 x I0~ 6 kg, m = 1 X 10 -8 fc#, F = 4s" 1 , 7 = 0.1s -1 , 
T = 295if and peaking frequencies shifted to the center 
of the plot. Curves from top to bottom are plotted from 
q = O.OliVm -1 , O.liVm -1 , liVm -1 , 10iVm -1 , lOOiVrn -1 . 



the pairs of peaking frequencies, one group clusters in 
the low frequencies and the other spreads out in the high 
frequencies. 

The low frequency group, the rescaled along peak cen- 
ter plot shown on the left of Fig[5l shares exactly the 
same peaking frequencies as those of the ancillary mass 
and we expect this behavior takes place when the target 
mass is resonating with the ancillary mass. This har- 
monic driven noise associates with the noise source F(uj) 
in Eq.fTrjJ and, as we argue before, does not contribute 
much to the overall thermal noise. The high frequency 
group spreads out to peak frequencies uj c = IQ.lkHz, 
10.5kHz, 14.2kHz, 33.3kHz and 101kHz for the varying 
spring constant g. The anharmonic noise with respect to 
the ancillary mass associates with the noise source f(ui) 
in Eq. (|14p . Our cooling mechanism in this identical case 
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FIG. 5: (Color online) Rescaled plot of DNS of the target 
resonator m with parameters: k — K — liVm -1 , M = 1 X 
W~ 6 kg, m = 1 X W~ S kg, F = 4s -1 , 7 = 0.1s -1 , T = 295K 
and peaking frequencies shifted to the center of the plot. The 
left and the right subplots corresponds to the low and the 
high peak frequency groups. Curves from top to bottom are 
plotted from g = O.OliVm -1 , O.liVm -1 , liVm -1 , lOiVm -1 , 
lOOiVrn -1 . 



is fairly effective, with the noise level reached down to 
4.03, 3.73, 2.72, 2.12 and 2.04 times a common factor 
of 10 -21 m 2 , respectively, identical to the values of the 
ancillary mass. 



B. The Large Detuning Case 

We assume, in this case, the spring constants K = 
lOOOiVm -1 and k = liVm -1 ; the masses M = 1 x 10 -3 fcg 
and m — 1 x 10 -8 fc(7; thus the natural oscillating fre- 
quencies for the two masses are retained. The damp- 
ing coefficients and the initial temperature are left un- 
changed. Varying the spring constant of the middle 
spring g over the same five values gives the noise spec- 
trum plot of the target mass m shown in FigJBJ The 
plot is again rescaled to the center along the peaking 
frequencies uo c = IQ.lkHz, 10.5kHz, 14.1kHz, 33.2kHz 
and 100.5kHz. 

We note that the noise peaks at one frequency for each 
value of the spring g. These peaking frequencies are close 
to those in the identical case above but the noise sup- 
pression rate, as we have expected, is much better. That 
means the large detuning not only helps suppress the 
noise source stemmed from the coupling mass M to neg- 
ligible amplitude but also makes the feedback more ef- 
fective for countering the target's noise. The 5 values of 
the spring constant g corresponds to noise fluctuations of 
4.03, 3.70, 2.03, 0.37 and 0.04 times a common factor of 
10 -21 m 2 . 

Fig[7| shows the plot of the ancillary resonator's noise 
spectrum, again rescaled to the center along the peaking 
frequencies. These peaking frequencies are very close to 
those of the target mass. We predict that the large de- 
tuning between the springs K and k puts the ancillary 
mass M into a particularly passive role that reflects the 
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V. NOISE IN MOMENTUM SPACE 
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FIG. 6: (Color online) Rescaled plot of DNS of the target 
resonator m with parameters: K = lOOOJVm -1 , k = lNm~ , 
M = 1 X W~ 3 kg, m = 1 X 1(T 8 %, T = 4s~\ 7 = 0.1s -1 , 
T = 295K and peaking frequencies shifted to the center of 
the plot. Curves from top to bottom are plotted from g = 
O.OlAm -1 , O.lTVm -1 , lTVm -1 , lOAm -1 , lOONm -1 . 



A. The Model with Harmonic Driving Force 

The effective temperature is usually denned according 
to the equipartition theorem through the equality be- 
tween the thermal energy and the mechanical potential 
energy. The intrinsic nature of our model renders the 
stored potential energy of the target and of the environ- 
ment in their shared link, the middle spring g, inseparable 
and thus forbids our discussion of the effective tempera- 
ture in the coordinate space. To make our model relevant 
to the discussion of cooling, we examine the system dy- 
namics in its momentum space and define the effective 
temperature according to the kinetic energy, i.e. the mo- 
mentum fluctuation, of the target. 

In addition, to follow the convention along previous 
literature, a harmonic driving force term 

G(t) = F e~^ ot , (27) 




co/Hz 

FIG. 7: (Color online) Rescaled plot of DNS of the ancillary 
resonator M with parameters: K = lOOOiVm -1 , k — \Nrn~ 1 , 
M = 1 x W~ 3 kg, m = 1 x W~ 8 kg, T = 4s _1 , 7 = 0.1s" 1 , T = 
295 AT and peaking frequencies shifted to the center of the plot. 
Curves from left to right are plotted from g = O.OlA^m -1 , 
CliVm" 1 , liVm -1 , lOiVm -1 , WONm' 1 . 



where loq is the driving frequency, should be added to 
the feedback loop (Cf. Figfl]). In our case, this harmonic 
force appears as an extra term in the equation of motion 
of the ancillary resonator M. 

Following the above arguments, we rewrite the coupled 
system of classical Langevin equations Eq.(fT]) and Eq.([2]) 
as 



p + jp + muj x — gQ = mf 



(28) 



P + TP + MVl 2 Q — gx = M(F + G(t)) (29) 

where p — mi and P = MQ denote, respectively, the 
momenta of the target resonator and the ancillary res- 
onator. The thermal fluctuations / and F still obey the 
same set of relations Eq. ([3]) and Eq. ([J| . 

Consequently, we have a modified linear response (Cf. 
Eq. (fl"4|) , Eq. (fl~5|) and Eq. (fT6)) ) after combining the system 
of equations in the frequency domain (see Appendix [D]) 

p(uj) — ium Lf(u))f(uj) + Lf{lo)[F(uj) 

+ 2ttF q S{lu + lu )}\ (30) 



minute motion of the target. This also helps explain why 
the noise suppression is especially effective in this case by 
the fact that the mass M has its speed comparable to the 
target m but in an opposite direction such that the force 
it exerts through the spring g can favorably counteract 
the motion of the target m. Nonetheless, the reinforced 
role that the ancillary mass M plays means that itself 
does not belong to the target system. As shown in the 
figure, the ancillary mass almost retains its original noise 
level throughout the varying values of the spring g. 



where Lf(uj) and Lp(uj) are the susceptibilities Eq.fT5]) 
and Eq. (fT6|) . We can thus arrive at the noise spectrum 
in momentum, after following the same routine of com- 
putations, 



S p (uj) 



1 

2^ 

= IknTm 



Auj' (p{u)p{uj')) 



(31) 



W 2 7 e 



(^ 



,2\2 



+ 



2nFQLolm 2 



S(uj + u> q )Lf(cj )L f (-u>o) (32) 
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B. The Effective Temperature 

The square of momentum of the target is defined sim- 
ilarly as in Eq. 



(P 2 ) 



1 

2^ 



du Sp(u>). 



(33) 



which can be separated into two parts: the variance of 
momentum, i.e. the fluctuation or noise in momentum 
space, 



(Ap) 2 \ = — - — m 



duj- 



W 2 7 e 



Wg — LO 2 ) 2 + LO 2 ^ 2 



(34) 



and the mean squared, i.e., the square of the steady state 
value of momentum, 

(p) 2 = FQU! 2 m 2 L F (uj Q )L F (-LLi ). (35) 

Note that the harmonic force term only contributes to 
the first moment of the momentum. 

Since thermal dissipation only induces variance of the 
target's momentum, the effective temperature of the tar- 
get can inversely be determine by the variance of the 
kinetic energy through the energy equipartition theorem 



T e s = 



(Ap)' 



(36) 



Eq. (|34p is independent of the two parameters Fq and ujq 
of the harmonic driving force. We conclude that the har- 
monic force bears no effect in reducing the system's effec- 
tive temperature and hence in cooling. Eq. (j34|) is also not 
analytically integrable, numerical integration shows the 
integral takes value similar to the case when the spring 
constant g — 0. That means the feedback is only effective 
in the displacement domain but not in the momentum 
domain. 



VI. DISCUSSION AND REMARKS 

In summary, we have proposed a theoretical model to 
demonstrate the general self-feedback type noise suppres- 
sion technique through the coupling between the target 
resonator and an adjuvant system. In particular, such 
a self-feedback is achieved through the dynamics of the 
adjuvant system. The explicit delay function has been 
given and its efficacy in reducing the noise spectrum ver- 
ified. We have also used numeral results to confirm our 
observations. 

Before concluding this paper, we add some remarks as 
follows: 

First, the theory of thermal equilibrium state can be 
used to explain why the system comprising two resonat- 
ing masses has its noise reduced to a limit when the 
spring constant g tends from zero to infinity. If g were 
zero, the stringing spring between the two masses would 



disappear and the system would comprise two indepen- 
dent resonators, each of which will attain its thermal 
equilibrium over time. In this case, the system has two 
degrees of freedom. According to the energy equipar- 
tition theorem, the system must contain twice as much 
energy as fcsT/2. For the other limit, when g tends to 
infinity, the two resonators become a rigid body and pos- 
sess only one degree of freedom, which means the system 
would contain only ksT/2 of energy. This explicit pic- 
ture clarifies how the combined system reduces its energy 
over the increasing value of g. 

However, how each resonator in the combined system 
reduces its own energy is an open question though we 
can intuitively think the heat bath may play crucial role 
in the asymmetric thermalization of the two resonators 
as tow open system. Another remark is that the energy 
equipartition theorem is applicable only after the system 
enters thermal equilibrium state. When interaction oc- 
curs between the two degrees of freedom contributed by 
the two masses, such an illustration is not appropriate. 
This is the reason why the nonlinear character of the 
noise suppression rate cannot be explained by the energy 
equipartition theorem alone. 

In addition, we have assumed there does not exist cor- 
relation between the two noise sources for the two masses, 
i.e. (f(t)F(t')) = 0. However, when the two masses oscil- 
late very closely with each other, the above assumption 
based on independent thermal environments will not hold 
and we need to consider the case where 



</ (W)> ¥> o. 



(37) 



Another setting we shall consider is when the motions 
of the two masses are quantized: we need to change the 
fluctuation-dissipation relations Eq.Q and Eq.((4]) to 



(/(t)/(tO) <* w(t,t',h), 

(F(t)F(t')) cx U(t,t',h). 



(38) 
(39) 



where W and U are not simply <5-functions and depends 
on the Planck constant H. We will present the general 
investigation on the case with quantum fluctuations and 
thermal bath correlations, but here we have concentrated 
on the simple case. 
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APPENDIX A: DERIVATION OF THE 
FEEDBACK DELAY FUNCTION 



Fourier transforming Eq.Q gives 



gx(u))/M + F{uj) 



Q 2 ~uj 2 



iluT 



(Al) 



where we use x, Q and F to denote the Fourier transforms 
of x, Q and F, respectively, 



x(u>) 



dtx{t)e- lu, \ 

) 

/CO 
dtQ(t)e~ iut , 
-OO 

dtF{t)e" luJt . 



(A2) 
(A3) 
(A4) 



integrate by parts with respect to (t — r) 

-2 



Q(t) 



V40 2 - r 2 Jo 



d(t - T)<f>(x,t) 



(A8) 



\nt-r) 



sm 



x exp 
1 



ft 2 v ' y V4fi 2 - r 2 

d</)(a;,r) 



< - r 



\/4ft 2 - r 2 



r 0(x,i) 
d(t - t) 



(A9) 



d(t - t) 



2 exp 



1 



r(t-r) 



exp 



±iV4fi 2 -r 2 (i - r) 



r - iV40 2 - r 2 



exp 



\Wm 2 - T 2 {t- T ) 



r + WAffl - r 2 

where the factor in the last two lines constitute the func- 
tion h{t — t) as defined in Eq. lfTTj) . 



Substituting the first two transform of the above into 
Eq. ljAlj) and taking the inverse Fourier transform of Q, 
we get 



2n J_ OQ (iV - u) 2 ) + iTlo 



The linear integral with respect to w can be viewed as a 
contour integral along a semicircl e in the up per half-plane 
with two poles at u>± = [iT ± \/4f2 2 — T 2 ]/2. Assuming 
r <C fi, the corresponding residues are, 



R±(t) = T 



1 



V4^ 2 - r 2 
i 



0(i- t) x 



(A6) 



exp 



where 9it — t) denotes the unit step function originated 
from the positive locus of the integration path. Using 
Cauchy's theorem, Q(t) is reduced to a single- integral 
form 



Q(t) = 



V40 2 - r 2 
i 



dr d>(x, t) x 



(A7) 



exp 



Tit r) 



-V40 2 -r 2 (t-r) 



where we have used the shorthand cj>{x,t) as defined in 
Eg. pop . To write in our desired form, we further 



APPENDIX B: DERIVATION OF THE 
RESPONSE FUNCTION AND THE NOISE 
SPECTRUM 

To find the response function, we first take the Fourier 
transform of Eq.Q and Eq.|T|) to get, respectively, 
Eq.(|AT]) and 



x(u>) 



gQjuj)/m + fjuj) 
uj 2 — uj 2 + iujj 



(Bl) 



where we let x, Q and / denote the Fourier transforms 
of x, Q and / as we did in Appendix [X] Substituting 
Eq.(|AT]) into Eq. fBTj) , we have 



x{ll>) 



g gx(w)/M + F(w) 
m fl 2 — uj 2 + iuiT 



+ f(u)j {Q 2 -u 2 +iun)~ 1 
(B2) 



Reshuffling the terms and putting again x (w) on one hand 
side, the equation becomes 



x(ui) 



fjuj) + gF (u)/[m(ft 2 -uj 2 + ijJT)\ 
Cu 2 — uj 2 + iuj^f — g 2 /[mM{Q 2 — uj 2 + iujT)] 



(B3) 

To find explicitly the response function, we group to- 
gether the real terms and the imaginary terms in the 
denominator to write i(o;) in a familiar form similar to 
the case where the coupling mass M were not present; 
whence the denominator becomes 
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O 2 



mM (ft 2 



,2\2 



-uj 2 r 2 



7 + 



mM in 2 - uj 2 ) 2 + lu 2 T 2 



(B4) 
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and we can define the terms in the two brackets as in 
Eq.CQJ and Eg. ([15]). 

Using Eq. (fT4|) . the associated fluctuation-dissipation 
relations of Eq. §5§ and Eq. (0J in the frequency domain 



/>)/>')) = 47rfc B T^<S(w+w / ), (B5) 
/ m 

F{lo)F(lo')\ = 4 w k B T—S{iJ + u'), (B6) 
/ ivi 



and the independence between the two noise sources, we 
find by using definition Eq. (|19j) 



S x {lo) = 2k B T 



7 r 

— L f(u)L f (~uj) + —L f (lo)L f (-lo) 



where Lf(u) and L f (lo) are defined in Eq. (|15p and 
Eq.dUl). Seeing that 



L F (u) 



g/m 



n 2 - lo 2 + iwY 



(B8) 



we can factor out Lf(Lo)Lf(—Lo)/m from the bracket in 
Eq.dEZD 



S x {u>) 



2k B T 



7 + 



g 2 T/mM 



(B9) 



1 



ft 2 - lo 2 + iuY fl 2 - uj 2 - iuoY 



2k B T 



1 



<J 



(lo 2 — lo 2 ) 2 + lo 2 ^/ 2 

2 r 



(BIO) 



mM (n 2 - lo 2 ) 2 + lo 2 Y 2 
which equals to Eq. (f2"0)) . 



APPENDIX C: DERIVATION OF THE 
VARIANCE OF DISPLACEMENT 

The integrand of Eq. (f2"2")) has highest order of lo 4 in the 
denominator and can not be integrated analytically. To 
make the integrand integrable, we introduce a parameter 
n(K) = g/(g + K) as a variable of the spring constant 
K and approximate the integrand by expanding Q 2 with 
respect to fi up to 2nd order in the vicinity of [i = 0, in 
other words, when K — > oo. That is, we let 



du| 

d// 



(Cl) 



To calculate the three terms in the expansion, we first 
substitute K with g(l — in Eq. (fTT)) . 



gM- 1 ^- 1 - lo 2 



mM (gM- 1 ^- 1 - lo 2 ) 2 + lo 2 Y 2 ' 



(C2) 



Since the lowest-order term in the numerator and the 
denominator is, respectively, /i" 1 and /x -2 , the second 



term of the above formula goes to when fi — » 0. Thus 
we find the first term in the expansion lo 2 |^_ = uj 2 . The 

first-order derivative of uj 2 (/i) reads 



d// 



-( 9 M 



-v- 1 



,2\2 



w 2 r 2 ] 



iM 2 



[(jM-V^-w^ + ^r 2 ] 1 



(C3) 



Again by comparing the coefficients of the lowest-order 
terms in the numerator and the denominator, we shall see 
the second term in the expansion diDg/d/i^^o = —g/ m - 
The third term in the expansion can be obtained similarly 
by observing the limiting behavior of the numerator and 
the denominator of the second-order derivative of uj 2 (/x) 

g 3 ^gM'^-^igM- 1 ^- 1 - lo 2 ) 2 



= lim 



ifi- 1 - lo 2 ) 2 +lo 2 T 2 Y 



dy 2 

(C4) 

which reads — 2(M/m)u> 2 after taking the limit. Com- 
bining the results, we have 



M 



-fj, 



1 2 
LO \1 . 



(C5) 



The truncation error rate introduced in this approxi- 
mation, by comparing Eq.( IC2[) and Eq. (|C5|l . is 



erf : 



Au, 2 



lo 2 - [lo 2 - (gfi - ML0 2 Li 2 )/m] 



-of 



(C6) 



Expanding lo 2 and [i gives 



erf = 



9 



Mlo 2 9 2 



g + K (g + K) 2 
x [{g + k)(F 2 + M 2 lo 2 T 2 ) - g 2 F] 



(F 2 + M 2 lo 2 T 2 ) - g 2 F 



(C7) 



where we have used the shorthand F = g + K — Mlo 2 . 
Since we only consider the usual cases with Y <C lo, the 
terms containing the damping coefficient can be omitted 



erf = 



[g 2 /(g + K) + Muo 2 g 2 /{g + K) 2 } F - g 2 
(g + k)F - g 2 



(C8) 



Expanding F and multiplying the numerator and denom- 
inator by (g + K) 2 , we find 

g 2 M 2 co 4 



{g + K) 2 [g(k + K - Mlo 2 ) + k(K - Mlo 2 



and at the other extreme of the expansion ji = 1 whence 
K = and g can take any value 



erf I 



M 2 lo 4 



K=0 



g(k - Mlo 2 ) - kMuo 2 ' 



(CIO) 



Therefore, the error can be sufficiently suppressed if we 
let g — ► oo and the expansion of Co 2 around /j, = 1 is 
validated. 

Following the same reasoning and procedure, we can 
approximate the effective damping coefficient Eq. (fT5f . 



7e 



7 + 



MY 



-M 



(Cll) 
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which becomes ^-independent. Substituting Eq. (|C5[) and 
Eq. (|Clip into Eq. (j2"2"|) , we get, after minor algebra, 



w 2 r i , P1 9 , 



e J — oo 



where 



, _ 7 + {MT/m)n 2 
7 ~~ {M/m)n 2 + 1 



oj 2 - {g/m)n 
{M/m)fi 2 + 1 



(C13) 
(C14) 



can be considered the effective damping coefficient and 
vibrating frequency of the mass m after the approxima- 
tion. We can apply a limiting process to Eq. (|C12j ) 



12 roo 



{x * (t)) = , iin 



-&0J 



miT"f e t^o J _ OQ [uj 12 — lo 2 ) 2 + o; 2 7' 2 

(C15) 

and compute the integral using the theorem of residues 
as we did in Appendix [X] 



(At)) 



12 



g + K 



g{k + K) + kK' 



(C16) 
(C17) 



APPENDIX D: DERIVATION OF EFFECTIVE 
TEMPERATURE 

Fourier transforming the Langevin equations in the 
momentum space Eq. (l28)) and Eq. ([29|) and recognizing 



p(u>) = iwx(ui) 
P(u>) = iu>Q(u>), 



(Dl) 
(D2) 



we have 



{lo 2 — lo 2 + iu^)f>(uj) — —P(cj) = iojmf(uj) (D3) 



M 



(n 2 -uj 2 + iujT)P(uj) - —P{lo) = iujM{F{uj) 



+ 2irF d(LO + w )). 

(D4) 



We can follow the lines in Appendix [Bl at this point to 
get the Fourier transform of the momentum of the target 
resonator 



mf(uj) + g 



F(uj) + 2wF S(uj + w ) 
51 — u> 2 + icoT 



oj 2 — uj 2 + ioj'-f - 



1 



mM n - uj 2 + iujT \ (U '' 1 

which can be written in a compact form as in Eq. ([30l) 
after grouping the terms. The density noise spectrum in 
momentum can then be computed 



S p (u>) 



1 

2tt. 
2kBTui 2 m 2 



duo' (p(u)p(cj')) (D6) 

r 

L f(u)L f (-Lj) + —L F (u)L F (-u) 

(D7) 



+ 2TrFQU! 2 m 2 S(uj + u!q)Lp(uj)Lf(—u)o) 



2kBTm- 



W 2 7e 



(w 2 — lu 2 ) 2 + UJ 2 "/ 2 
+ 2nF 2 uj 2 m 2 5(uj + oj )L f (uj )L f (—uj ) 



(D8) 
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